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^ ' Abstract 
>v 

The well known DuhamePs principle allows to reduce the Cauchy problem for a linear in- 
homogeneous partial differential equation to the Cauchy problem for the corresponding homo- 
geneous equation. In the paper one of the possible generalizations of the classical Duhamel's 
principle to the time-fractional pseudo-differential equations is established. 

: 1 Introduction 

The role of the classical "Duhamel's principle", introduced by Jean-Marie-Constant Duhamel in 
1830th, is well known. The main idea of this famous principle is to reduce the Cauchy problem 
for a given linear inhomogeneous partial differential equation to the Cauchy problem for the cor- 
responding homogeneous equation, which is more simpler to handle. In this paper we establish 
a fractional analog of the Duhamel's principle with respect to the following Cauchy problem for 
inhomogeneous time-fractional pseudo-differential equations 

' D«u(t,x) = A{D x )u(t,x) + f(t,x), t>0, x£R n , (1) 

^(0, x) = <p k {x), x G R n , k = 0, m - 1, (2) 

where a E (m — m > 1 is an integer; f(t,x) and (fk{x), k = 0, m — 1, are given functions 

in certain spaces defined later; D x = (D±, D n ) Dj = — j = l,...,n; A(D X ) is a 

pseudo-differential operator with a symbol A(£) defined in an open domain G C R n ; and is 
the operator of fractional differentiation of order a in the Caputo sense (see, for example, 0111121) 

r fit), i a = o, 

(#?)/(*)= ITlMo0F> 0<a<1 ' (3) 
{ D T m (fif{v)), m<a<m+l, m=l,2,.... 

Note that at the same time other generalizations of the Duhamel's principle based on other defini- 
tions of fractional derivatives are possible, but in this paper we do not stop by these possibilities. 

In our analysis we essentially relay on the results obtained in the paper [Hj, where the Cauchy 
problem for fractional order homogeneous pseudo-differential equations of arbitrary order a (a > 0) 
is studied. Note that fractional order inhomogeneous equations can not be directly reduced to 
the corresponding homogeneous equations and the classical Duhamel's principle is not applicable. 
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To our best knowledge, solution of the Cauchy problem for general inhomogeneous fractional or- 
der pseudo-differential equations require construction and estimation of the corresponding Green's 
function or combination of classic Duhamels principle and some integral equations jU |5J El El © • 
The fractional Duhamel's principle established in the current paper can be applied directly to in- 
homogeneous fractional order equations reducing them to corresponding homogeneous equations, 
at least, in the framework of the Cauchy problem lfl |) . 

This paper is organized as follows. In Section 2 we recall some facts related to the fractional Ca- 
puto derivative, the classic Duhamel's principle, and the pseudo-differential operators with singular 
symbols. In Section 3 we formulate the main result of the present paper, namely the fractional 
analog of the Duhamel's principle. In this section we also demonstrate some applications of the 
obtained principle. 



2 Preliminaries 

2.1 On the Caputo fractional derivative 

It follows from the definition @ of the Caputo fractional derivative of order a € (m — l,m] for a 
given function f(t), that its m-th derivative f^ m \t) has to exist and satisfy certain integrability 
conditions. We require throughout the paper that these conditions are fulfilled. We note also that 
if a = m, where m > 1 is an integer number, then D°f(t) = f m (t), which immediately follows 
from the definition. 

Further, denote by J 7 , 7 > 0, the fractional order integration operator 

J1/(t) = fMi ,(t - T)7 " I/(T)<iT ' 

with J° = I, I is the identity operator. It is well known [2], that the family {J 7 , 7 > 0} possesses the 
semi-group property. Namely, J^ 1 ^ 2 = JT 1 J^ 2 = J 72 J 71 , j t > 0, 72 > 0. The Caputo derivative 
can be written in the form (see 

D»f(t) = J m - a /( m )(t), m - 1< a < m. (4) 

If a function f(t) satisfies the condition 

/( fc )(0) = 0, k = Q,...,m-l, (5) 

then for a £ (m — 1, m] and a real (5 > we have 

jP+ a f(t) = jP +m D™- a f(t). (6) 

The relation © is obviously fulfilled, if a = m. Assume m — 1 < a < m. Then for f(t), satisfying 
© the equality f(t) = J m f {m) {t) holds. Hence, 

jP +a f(t) = J^ +a J m f^ m \t) = J p+m J a f( m \t). (7) 

From the other hand D™~ a f(t) = J™-{™-a) /M ( t ) = This fact and the re i a ti on (J7J) 

imply (JBJ). 

Further, if / does not satisfy the conditions ©, then still the relation analogous to © holds, but 
with D™~ a , the fractional derivative in the sense of Riemann-Liouville, instead of D™~ a . Namely, 
in this case the relation 

jP +a f(t) = jP +m D™~ a f(t), (8) 
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holds, where 

1 d m r l f(r) 

D\f(t) = r^— / -, m - 1< 7 < m, 

+JKJ r(m-7) dt m Jo (t - r)T +1 - m dr 

and = / m (t). Indeed, in this case we have f(t) = J m f {m) (t) + J2T=o j£ tP**> which y ields 



j(3+m ja Am) / j.\ , /( fe )(Q)t fc +£+" 



k=0 

From the other hand using the relation (see |2j) 

D a + f{t) = D«f(t) + J2 r ,/ 1^ ' m - 1< a < m, 

r ( fc " « + 1) 

we have 

2.2 Basic spaces of functions and pseudo-differential operators 

In this section we briefly recall some necessary notions and facts, which we use in this paper 
referring the reader to [31 EH; f° r details. Let G be an open domain in R n and the system of open 
sets {gk}kLo be a locally finite covering of G, i.e. 

oo 

G = (J g k , 9k CC G. 

k=0 

Let any compact set K C G has a nonempty intersection with finitely many sets g k . Denote by 
{^fclfcLo a sm °°th partition of unity of G. Further, let 1 < p < oo and a function f(x) be in L p (R n ) 
whose Fourier transform F f has a compact support in G. The set of all such functions endowed 
with the convergence defined in Definition 1 is denoted by ^c,p{R n )- 

Definition l.A sequence of functions f m G ^Q p [R n ), m= 1, 2, 3, ... is said to converge to an 
element f £ ^ G , P ( Rn ) i S : 

1 ) there exists a compact set K C G such that the support supp Ff m C K for all m € N; 

i 

2) the norm \\f m - f \\ Lp = (f Rn |/ m (z) - f {x)\ p dx)p -> for m -> oo. 

According to the Paley- Wiener-Schwartz theorem elements of ^G jP (i? n ) are entire functions of 
exponential type, which, restricted to R n , are in the space L p (R n ). 

The space ^fa,p(R n ) can be represented as an inductive limit of some Banach spaces. Namely, 

let 

N N 

G N =\Jg k , x*(0 = EW 

k=l k=l 

Denote by *$>n the set of functions / £ L p (R n ) satisfying the following conditions: 

a) suppFf c Gat; 

b) suppFf fl supp 6j = for j > N; 
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c)p N (f) = \\F- 1 X N Ff\\ Lp <™. 
Here by F 1 we denote the operator inverse to the Fourier transform F. It is not hard to verify 
that (see 0113) 

$G„ = md lim 

Let A(D) be a pseudo-differential operator with a symbol -A(£), which is a real-analytic 1 function 
in G. Outside of G or on its boundary A(£) may have singularities of arbitrary type. For a function 
<p(x) € ^ , G,p(-R n ) the operator A{D) is defined by the formula 

As is shown in ^3 EH the space ^Q >p (R n ) is invariant with respect to the action of such pseudo- 
differential operators and these operators act continuously. 



2.3 The classic Duhamel's principle 

Recall the classic Duhamel's integral and Duhamel's principle. The Duhamel's integral (see, e.g. 
|13[ I14j ) is used for representation of a solution of the Cauchy problem for a given inhomoge- 
neous linear partial differential equation with homogeneous initial conditions via the solution of the 
Cauchy problem for the corresponding homogeneous equation. Consider the Cauchy problem for 
the second order inhomogeneous differential equation 

-^(t,x)=Lu(t,x) + f(t,x), t>0,xeR n , (9) 
with homogeneous initial conditions 

u(0,x)=0, ^{0,x)=0, (10) 

where L is a linear differential operator with coefficients not depending on t, and containing 
the temporal derivatives of order, not higher then 1. Further, let a sufficiently smooth function 
v(t, t,x), t>T, t > 0, x E R n , be for t > r a solution of the homogeneous equation 



-Q^(t,T,X) = Lv(t,T,x), 



satisfying the following conditions: 



dv 

v(t,T,X)\t=r = 0, —(t,T,x)\ t=T = f(T,x). 

Then a solution of the Cauchy problem Q, (|10j) is given by means of the Duhamel's integral 

u(t,x) = / v(t,T,x)dr. 
Jo 

The formulated statement has been entitled as the "Duhamel's principle". 

1 This condition can be essentially weakened. See the construction in jlll I12| . where continuous symbols are 
considered. 
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An analogous construction is possible in the case of the Cauchy problem with a homogeneous 
initial condition for the first order inhomogeneous partial differential equation 

— (t,x) = Mu(t,x)+f(t,x), t>Q,xER n , 

where M is a linear differential operator containing only spatial derivatives, and with coefficients 
not depending on t. 



2.4 The Duhamel's principle for integer a = m, m > 1 

Consider the Cauchy problem (1), (2) in the case of integer a = m > 1, = J^-, i.e. 

(t, x) = A(D x )u(t, x) + f(t, x), t > 0, x E R n , 
(0,x) = <fk(x), x E R n , k = 0, m — 1. 



dt 

d k u 



dt k 

In this case the Duhamel's principle is formulated as follows. Let U(t, r, x) be a solution of the 
Cauchy problem for a homogeneous equation 

^ = A(D X )U, 0<r<t, (11) 

d k U 



dt k 



(t,T,x)\ t=T = 0, k = 0,...,m-2, (12) 
^r(t,r,x)\ t=T = f(r,x). (13) 

Then the function 

u(t,x) = / U(t,T,x)dr (14) 
J 

is a solution of the Cauchy problem 

A(D x )u = f(t,x), (15) 



k 



dt r 
u 

(0,s) =0, k = 0,..., m-1. (16) 



dt k 

The proof of this statement can be found, for instance, in [7]. However, for the completeness, 
we reproduce the proof. 

Obviously u(0, x) = 0. Further, for the first order derivative 

du f f d 

— (t,x) = U(t,t,x)+ —U(t,T,x)dT, 

it follows from (JT2J) that ^f(0,x) = 0. Analogously we calculate |jf(0, x) = 0, k = 1, ...,m - 2. 
Consequently, for the derivative of (rn — l)-th order 

Qm—\ Qm—2 j-t Qm—1 

— T u(t,x) = — ~U(t, t, x) + / T U(t,T,x)dr, 

g t m-i v ' > dt m ~ 2 v ; Jo dt™- 1 v ; 
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we obtain g f m-i u (0, x) = 0. Therefore, the function u(t, x) in (|14|) satisfies the initial conditions 
Q16JI . Moreover, 



d m u , . 3 



a" 1 r r 

A{D x )u = — J U(t,T,x)dT-A(D x ) J U(t,T,x)dT~- 

T U(t,t,x) + J —U(t,r,x)dT - J A(D x )U(t,T, X )dT 



dt 

rt Qrn 



= /(*, x) + J^ [g^U(t, t, x) - A(D x )U(t, r, x)\dr = f(t, x). 
Hence, u(t, x) in (j!4|) satisfies the equation (|T5)l as well. 



2.5 The representation formula for a solution of the Cauchy problem for ho- 
mogeneous fractional order equations 

Now we consider the Cauchy problem ((H). (J2J) for arbitrary a > 0. Note that in this case Duhamel's 
principle can not be applied directly. For the Cauchy problem ({! ]) . 1(2 )1 in the homogeneous case 
(i.e. f(t,x) = in Equation (fT|)) the following representation formula for a solution was obtained 
in 0: 

=53^~ 1 ^a(* a A(Z) ai )) Vfc _ 1 (a!), (17) 
fc=i 

where J fc is the fc-th order integral operator, E a (t a A(D x )) is a pseudo-differential operator with 
the symbol E a (t a A(£)) and i£ (2) is the Mittag-Leflier function (see [T5| ) 



3 Main results 

3.1 A fractional Duhamel's principle in the case < a < 1 

Assume < a < 1. First we formulate a formal fractional analog of the Duhamel's principle and 
then we show how to apply this principle in various situations. 

Theorem 1. Suppose that V(t,r,x), < r < t, x € R n , is a solution of the Cauchy problem 
for homogeneous equation 

D%V(t, r, x) - A(D x )V(t, r, x) = 0, t > t, x€ R n , (18) 

V(t,t,x) =Dl- a f(r,x), x€R n , (19) 
where f(t,x) is a given function satisfying the condition f(0,x) = 0. Then 

v(t,x) = f V(t,T,x)dr (20) 
Jo 

is a solution of the inhomogeneous Cauchy problem 

D«v(t, x) - A(D x )v(t, x) = f(t, x), (21) 
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v(0,x) = 0. 



(22) 



Proof. Notice that in accordance with (|T7)l a solution of the Cauchy problem qi8 |) .()19| ) is 

represented in the form 

V(t,T,x)=E a ((t-T) a A(D x ))Dl- a f(T,x). (23) 

Further, apply the operator J a to both sides of Equation (|2Tj) and use the relation J a D^v(t, x) = 
v(t, x) — v(0, x), to obtain 

v(t,x) - J a A{D x )v{t,x) = J a f{t,x). 
A solution of the last equation can be represented as 

oo 

v(t,x) = J2J an+a A n (D x )f(t,x). 

n=0 

It follows from (JSJ) (with /3 = an and m = 1) that for arbitrary function g(t) satisfying the condition 
g(0) = 0, there holds J an+a g{t) = J an+1 D\~ a g(t) . Taking this into account we have 

t 00 



,(t,x) = £ J^A"(D x )Dl-f(t,x) = f £ {t r T ( rn ^ x) Dl-f(t,x 

n=0 J0 n=0 1 {an + l > 



f E a {{t - T ) a A(D x ))Dl- a f(t, x). (24) 
Jo 

Comparing (P3J) and (PH we arrive at lj20)l. 
Remark. 

1. If a = 1, then Theorem 1 coincides with the classic DuhameVs principle, because in this case 
Dl- a f(t,x) = f(t,x). 

2. It is well-known that the fractional derivative D\~ a f{t), < a < 1 exists, if f(t) G ^4C[0 < 
t < T], where T is a positive finite number and AC[0,T] is the class of absolutely continuous 
functions (see 



3. The condition f(0,x) = in Theorem 1 is not essentially restrictive. For arbitrary f(t,x) in 
the formulation of Theorem the Cauchy condition Uty) has to be replaced by 



V(r,T,x) = D^ a f(T,x), xGiT, 

where D^T a is the operator of fractional differentiation of order 1—a in the Riemann-Liouville 
sense. 



Denote by C (m) [t > 0; ^ Gt2 (R n )] and by AC[t > 0; ^ Gt2 (R n )] the space of m times contin- 
uously differentiable functions and the space of absolutely continuous functions on (0; +oo) with 
values ranging in the space ^c,2(R n ), respectively. 

Theorem 2. Let ip (x) G ^ G , 2 (i? n ), f(t,x) G AC[t > 0; *G, 2 (# n )], Dl~ a f{t,x) G C[t > 
0; X I , G,2(-R n )] and f(0,x) = 0. Then the Cauchy problem (1), (2) (with < a < 1) has a unique 
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solution u(t,x) G C^[t > 0; ^ G2 (R n )] n C i f ^ °5 *G,2(^ n )]- ™s solution is given by the repre- 
sentation ^ 

u(t,x) = E a (t a A(D x )) Vo (x) + f E a ((t-T) a A(D x ))Dl~ a f(T,x)dT. (25) 

Jo 

Proof. The representation ()25j) is a simple implication of ()17|) and Theorem 1. The first term 
in (|23|) is studied in [2j in detail. Denote by v(t,x) the second term in (|23|) . For a fixed t > 
making use of the semi- norm of ^ at we have 

p 2 N (v(t,x)) = WF-^sM'L = \\x N Fv\\l 2 = 
= I \XA0\ 2 ■ I f E a ((t-rrA(0)FDl- a f(r,Odr\ 2 dt 

jR n Jo 

For Xjv(£) there exists a compact set K N C G such that suppx N (Q C K N . By using Cauchy- 
Bunjakowski's inequality we get the estimate 

p 2 N (v(t,x)) < [ \xA0\ 2 - f\E a {(t-TYA(m 2 dT- flFDl-y^Ofdrdt. 
Jk n Jo Jo 

The function f \E a ((t — r) a A{£,))\ 2 dr is bounded on K N . Consequently, there exists a constant 
C N > such that 

p 2 N (v(t,x))<C N [ \x N (0\ 2 - f\FDl- a f(r } i)\ 2 dTdi< 
Jk n Jo 

<C N f I \X N {0? ■ \FDl- a f(r,0\ 2 dtidT = 
Jo JR n 

= C N f a Wx N {t)FD l - a f(T } 0\\ 2 L 2 dT = C N j\ 2 N (Dl- a f(r,x))dr. 

It follows from the condition Dl- a f(t,x) G C[t > 0; ^ G , 2 (R n )} that the function p N ( D l~~ a f( T , x )) 
is continuous with respect to r G (0; t) and for a fixed t > and some N% the estimate 

p 2 N (v(t,x)) < C N ■ t • sup p 2 N (Dl- a f(r,x)) <C Ni -t- sup p 2 (f(r,x)) 

holds. Consequently, for every fixed t G (0; +oo) the function v(t, x) in (|2U|) belongs to the space 
^G,2{R n )- The analogous estimate is valid for §- t v{t,x). Thus v(t,x) G C^[t > 0; ^ Gt2 {R n )] n 
C[t > 0; ^ Gi2 (R n )]. Hence, v(t,x) G C^[t > 0; ^ G ,2(# n )] n C[t > 0; ^ G:2 (R n )], as well. The 
uniqueness of a solution follows from the representation formula for a solution of the homogeneous 
Cauchy problem. 

3.2 A fractional Duhamel's principle in the case of arbitrary a > 

Now we consider the Cauchy problem (1), (2) for arbitrary order a, m — 1 < a < m £ N. 
Obviously, in this case < m — a < 1. 

Theorem 3. Assume m > 1, m — 1 < a < m, and V(t,T,x) is a solution of the Cauchy 
problem for the homogeneous equation klty with the Cauchy conditions 

d k V 

— j -(t, T ,x)\ t=T = Q, k = 0,...,m-2, (26) 
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^-^(t,T,x)\ t=T = Dr a f(r,x), (27) 

where f(t,x), t > 0, x G R n , is a given function satisfying the conditions 9 -^r^ = 0, fc = 0, ...,m — 
1 . T/ien ^ 

v(t,x) = { V(t,r,x)dT (HFJ 
Jo 

is a solution of the Cauchy problem for the inhomogeneous equation \21\l with the homogeneous 
Cauchy conditions 

_(0,a;)=0, fc = 0,..,m-l. (28) 

Proof. It follows from the representation formula (|17j) that 

V(t, r, x) = J m ~ l E a {(t - r) a A(D x ))DT~ a f(r, x) (29) 

solves the Cauchy problem for Eq. (|18|) with the initial conditions (|26|) . (|27|). 
Further, apply the operator J a to both sides of the equation ()21l) and obtain 

m— 1 .j j rr\ \ 

v(t,x)-J2 ~ V -J a A(D x )v(t,x) = J a f(t,x). (30) 
i=o ^ 

Taking into account the conditions (|28|) . we rewrite Eq. ()3(J|) in the form 

v(t,x) - J a A(D x )v(t,x) = J a f(t,x). 
A solution of this equation is represented as 

oo 

v(t,x) = Y,J an+a A n (D x )f(t,x). 

n=0 

It follows from Q (with (5 = an) that for arbitrary function g(t) satisfying the conditions <?(0) = 
0, g'(fi) = 0, gl™- 1 ) (0) = 0, there holds J an+a g(t) = J an+l J m ~ l (D™- a g{t)). Taking this into 
account, we have 



v(t,x) = ]T J an+l J m ~ l A n {D x )DT- a f{t,. 



n=0 



Jo f-i. T(an + 1) 



n=0 

r-f 



/ r-^^^-r)^^))^/^^)- (31) 
Jo 

Comparing 1|29|) and (|31|) we obtain (J20J), and hence, the proof of the theorem. 

Remark. Note that if a = m, then Theorem 3 recovers the known classic Duhamel's principle 
we mentioned above (see, subsection 2.4)- Moreover, the conditions 9 ofk'^ = 0, k = 0, ...,m — 1, 
which we required in Theorem, are not essential. For arbitrary f(t,x), as a consequence of the 
relationship the formulation of the fractional Duhamel's principle takes the following form. 
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Theorem 4. Assume m > 1, m — 1 < a < m, and V(t,r,x) is a solution of the Cauchy 
problem for the homogeneous equation M8\) with the Cauchy conditions 

^(t,T,x)\ t=T = 0, k = 0,...,m-2, (32) 

£\7Th — 1 T T 

- w - T (t,T,x)\ t=T = D™- a f{T,x), (33) 

where f(t,x), t > 0, x G R n , is a given function. Then v(t,x) defined in 120\) is a solution of the 
following Cauchy problem for the inhomogeneous equation 

D?v(t,x)-A(D x )v(t,x) = f(t,x), 
d k v 

_(0,a;)=0, fc = 0,..,m-l. 



Theorem 3 and Theorem 4 allow to generalize the result of the paper ^ for arbitrary a > 0. 

Definition 2. A function u(t,x) G C {m \t > 0; ^ Gt2 (R n )] n C( m-1 )[t > 0; ^ Gt2 (R n )\ is called 
a solution of the problem (1), (2) if it satisfies the equation (1) and the initial conditions (2) 
pointwise. 

Theorem 5. Let ip k {x) G ^ Gi2 (R n ), k = 0,..., m - 1, f(t,x) G AC[t > 0; ^ Gi2 (i? n )], 

D™- a f(t,x) G C[t > 0; #G, 2 (-R")] and = 0, fc = 0,...,m - 1. T/ien the Cauchy problem 

(1), (2) has a unique solution. This solution is given by the representation 

m 

u(t,x) = Y, J k ' 1 E a (t a A(D x ))<p k _ 1 (x)+ 

k=l 

f J m - 1 E a ((t-rrA(D x ))DT~ a f(r,x)dr. (34) 
Jo 



Proof. Splitting the Cauchy problem (|T ]l .([2 ^1 into the Cauchy problem for the equation (^Q) 
with the homogeneous initial conditions and the Cauchy problem for the homogeneous equation 
corresponding to (0) with the initial conditions ((21), and applying Theorem 3 and representation 
formula (|17|). we obtain (|34j). The fact that 

m 

YJ^EUeAiD^.^x) G C^[t > 0;^ G , 2 (R n )]nC^[t > 0; * Gj2 (i? n )] 

k=l 

is proved in .3; . Further, since the m — 1-th derivative with respect to t of the last term in (|34[) 
belongs to j4C[[0, T]; ^ Gt2 (R n )] 2 , then the estimation obtained in the proof of Theorem 2 holds in 
this case as well. 

Remark. If f(t,x) and its derivatives up to order m — 1 with respect to t do not vanish at 
t = 0, then in accordance with Theorem 4 the representation formula \34\) takes the form 

m 

u{t,x) = J2 J k - 1 E a (t a A(D x ))<p k _ 1 (x)+ 

k=l 

2 T is an arbitrary positive finite number. 
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/* J m - l E a ({t - r) a A(D X )) D™~ a f (t, x)dr. 
Jo 

3.3 Examples 

1. Let < a < 1 and f(t,x) be a given function, f(0,x) = 0. Consider the Cauchy problem 

D?u(t, x) = Au(t, x) + f(t, x),t>0,xe R n , 
u(0,x) = ipo(x). 

In accordance with the fractional Duhamel's principle (Theorem 1) the influence of the external 
force f(t, x) to the output can be count as 

D?V(T,t,x) = AV(r,t,x), t>r,xeR n , 

V(r,r,x)=Dl- a f(r,x). 
The function V(t, t, x) = E a ((t — r) Q A)D^ _a /(r, x) solves this problem. Hence, 

u(t,x) = E a (t a A)<p (x) + f E a ((t-T) a A)Dl- a f(T,x)dT. 

Jo 

2. Similarly, if 1 < a < 2, and F(t, x) is a given function, which describes the outer force, then 
we deal with the Cauchy problem 

D°u(t,x) = Au(t,x) + F(t,x), t>0, xeR n , 

u(0,x) = (fo(x), u t {0,x) = ipiix). 

Again in accordance with the fractional Duhamel's principle (Theorem 4) the influence of the 
external force F(t, x) to the output can be count as 

D°V(r,t,x) = AV(r,t,x), t>r,xeR n , 

V(t,t,x) = 0, -q£(t,t,x) = D 2 + - a F( T ,x). 
The function V(r,t,x) = JE a {(t - T) a A)D 2 + a F{T,x) solves this problem. Hence, 

u(t,x) = E a (t a A)<p (x) + JE a {t a A)i Pl {x)+ 

f JE a ({t - T) a A)D 2 + - a F(r,x)dT. 
Jo 
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